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SUMMARY 


A limited experimental investigation has been conducted to confirm the 
existence of exhaust backflow from a cluster of three widely spaced nozzles 
operating in a near-space environment and to determine pressures and heat- 
transfer coefficients in the region washed by the backflow. Experiments were 
conducted in a 6l- foot-diameter vacuum sphere using a single solid-propellant 
rocket motor and a reflection plate. The reflection plate was used to simulate 
the Interference planes that would normally exist for a cluster of three noz- 
zles. Data were obtained for nozzle spacings of 8.02 and 17*05 nozzle exit 
diameters for pressures which simulated a nominal altitude range of 5.2 X 10^ 
to 1.8 x 10^ feet. The pressure data indicated that the backflow was super- 
sonic for several nozzle exit diameters ahead (upstream) of the nozzle exit 
plane and that the flow over small-diameter thermocouples in this region was 
in the free-molecule region. Average heat-transfer coefficients obtained for 
small-diameter cylinders (thermocouples) in this region were on the order of 
1 x 10”^ Btu/sq ft-sec-°R. 


INTRODUCTION 


The grouping of several nozzles in a cluster on a rocket-powered vehicle 
introduces a flow phenomenon that is a direct result of exhaust plume expansion. 
This flow phenomenon is generally referred to as exhaust backflow and occurs 
when the expanding exhaust plumes from adjacent nozzles impinge upon each other 
and create a shock pattern and attendant pressure gradient that cannot be nego- 
tiated by a portion of the flow in the plume boundaries. Consequently, this 
portion of the exhaust flow turns and flows in a direction opposite to the main 
nozzle flow. The high stagnation temperature of this backflow imposes a con- 
vective heat load on the vehicle base and nearby components and, unless suitable 
protection is provided for these items, destructive damage to the vehicle may 
result. Therefore, a number of investigations have been conducted (refs. 1 to 6, 
for example) to determine base heating rates for various configurations over a 
range of altitudes. Much of this work was experimental, as the many factors 
involved make analytical study very difficult, and is characterized by nozzle 


spacings less than the nozzle exit diameter and ambient pressures which simulate 
a maximum altitude of about 1.5 X 10^ feet. 

In space, where the exhaust plume expansion is a maximum, interference 
between exhaust plumes from a cluster of nozzles may occur even with large sep- 
aration between adjacent nozzles. The open-framework type of construction used 
for certain types of space vehicles would expose certain components to exhaust 
backflow from clusters of propulsive or control rocket motors. An example of 
this is the Surveyor spacecraft where the three vernier control rockets have a 
spacing of about 11.8 nozzle exit diameters and a base, in the normal sense, is 
absent. Even with this large nozzle spacing the exhaust plumes will impinge 
upon each other only a short distance from the nozzles and the expected backflow 
would wash components and structure located behind the control rockets. Although 
the heat transferred to these items from any backflowing exhaust gases should be 
small because of the low density of the backflow, it may become larger than 
desired if the backflow velocity is high or if backflow chokes escape passages 
and thus creates regions of relatively high pressure. In addition, the backflow 
may tarnish highly polished surfaces and thus increase the solar radiative heat 
flux to these surfaces. 

The exploratory investigation reported herein was undertaken to confirm the 
existence of exhaust backflow for large nozzle spacings in a near-space environ- 
ment and to investigate the pressures and heat transfer in the region washed by 
the unrestricted backflow. Tests were conducted in a 6l- foot-diameter vacuum 
sphere using a single solid-propellant rocket motor and a reflection plate. The 
reflection plate was used to simulate the interference planes that would normally 
exist for a cluster of three nozzles. Pressures and heat-transfer coefficients 
in the backflow region were determined for nozzle spacing ratios (ratio of dis- 
tance between nozzle exit diameters to the nozzle exit diameter) of 8.02 and 
17.05 over a nominal altitude range of 5*2 X 10^ to 1.8 x 10^ feet. 
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SYMBOLS 

area, for thermocouples and for calorimeters, sq. ft 

specific heat, Btu/lb-°F 

shortest distance between inside edges of adjacent nozzles measured 
in plane of nozzle exits, in. (see fig. 1) 

nozzle spacing ratio 

diameter of thermocouple wire or diameter of calorimeter disk, in. 

q j 

convective heat-transfer coefficient, - , Btu/sq ft-sec-°F 

Tt - T 

length of thermocouple wire loop, ft 
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backflow Mach number 

Nusselt number with thermal conductivity determined at stagnation 
temperature of backflow 

pressure, lb/sq in. absolute 

nozzle total-pressure ratio 

reflection-plate static-pressure ratio 

backflow total-pressure ratio 

heating rate, Btu/sec 
heat flux, Btu/sq. ft- sec 

gas constant, ft/°F 

backflow unit Reynolds number, l/ft 

Reynolds number used in correlating heat-transfer data, based on 
diameter of thermocouple wire with viscosity of backflow evaluated 
at stagnation conditions and density at free-stream conditions 

thermocouple Reynolds number based on wire diameter 

time, sec 
temperature, °F 


Tl, T2, etc. thermocouple or calorimeter position designation on reflection 
plate (see fig. 3(b)) 

w specific weight, lb/cu ft 

W weight, lb 

y distance perpendicular to surface of reflection plate, in. 

a angle between tangent to exhaust boundary at nozzle exit and nozzle 

center line, deg 

J3 angle between tangent to exhaust boundary and reflection plate at 

the point of intersection of exhaust boundary and reflection 
plate, deg 
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ratio of specific heats 


Subscripts: 

c combustion chamber 

e nozzle exit 

C heat meter or calorimeter 

s reflection-plate surface 

t total 

T thermocouple 

oo ambient or sphere conditions 

APPARATUS AND METHODS 


Facility 

The investigation reported herein was conducted in the Langley 61-foot- 
diameter vacuum sphere. The sphere pressure before firing of the rocket motor 

was approximately 3 X 10 millimeters of mercury for all tests and the post- 
fire pressure was approximately 0.26 millimeter of mercury. The altitude range 

corresponding to these pressures was about 3*2 X 10^ to 1.8 x 10^ feet. Initial 
evacuation of the sphere to a few millimeters of mercury was accomplished by 
positive displacement type of pumps and final pump-down by diffusion pumps. A 
McLeod gage was used to determine the sphere pressure before firing of the 
rocket motor. Post-fire sphere pressure was measured by an electronic vacuum 
gage. 


Model 

To eliminate the problem of obtaining equal performance and simultaneous 
ignition of three rocket motors and to reduce the sphere pressure rise to a 
minimum, a reflection plate was placed adjacent to a single rocket motor to 
simulate the interference planes that would normally exist for a cluster of 
three nozzles. This procedure is illustrated in figure 1. The l/4-inch- thick 
steel reflection plate was supported 10 feet above the floor in the sphere by 
a steel framework (fig. 2). The rocket motor and its combination holder and 
thrust measurement transducer (load cell) were supported adjacent to the plate 
by a section of steel pipe. The center of the nozzle was located 6 inches 

= 8.02^ and 12 inches = I 7 .O 5 ) from the plate being measured as shown 

in figure 2. 
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The polyurethane propellant rocket motor had a nominal thrust of 100 pounds 
and a useful burning time of about k seconds- The combustion-chamber pressure 
and temperature were approximately 650 lb/sq in. abs and b800 °F, respectively. 
The conical nozzle was an integral part of the motor casing and had a half-angle 
expansion of 30 with a nominal throat diameter of 0.3 inch and an exit diam- 
eter of 1.33 inches. For these dimensions and a ratio of specific heats of 
1.22, the exit Mach number was 3*83* A conical nozzle was used so that the 
angle through which the exhaust flow turned at the nozzle exit for the condi- 
tions of this experiment would more nearly approximate that which would occur in 
space for a contoured nozzle. (See ref. 7- ) 


Instrumentation 

The location and designation of the 1 6 thermocouple or calorimeter 
mounting holes and the six 0 . 060 -inch-diameter static-pressure orifices on the 
reflection plate are shown in the drawings presented as figure 3* Since the 
height of the thermocouples and calorimeters above the reflection plate and the 
location of the total-pressure probe was varied for each nozzle spacing ratio, 
configuration numbers were assigned to each motor firing. These configuration 
numbers and the corresponding thermocouple and calorimeter heights and probe 
locations are listed in the table of figure 3(*>)- The manner in which the ther- 
mocouples and calorimeters were installed in the reflection plate and the height 
designated is shown in the sketch of figure 3(a). 

Pressures . - In addition to the six static pressures measured on the plate, 
the sphere pressure, rocket-motor combustion- chamber pressure, and the total 
pressure of the exhaust backflow were also recorded. The plate static pressures 
and the backflow total pressure were measured with 0.1 lb/sq in. differential, 
inductance -type transducers calibrated for a pressure range equivalent to 
2 inches of water. These transducers were connected to the pressure orifices 
by 2-inch- long, l/8-inch inside diameter tubing. The sphere pressure was meas- 
ured with a similar type of transducer calibrated for a pressure range equiva- 
lent to 1 inch of water. This transducer was referenced to a small pressure 
chamber that was maintained at the pressure existing in the sphere prior to 
motor ignition. All other pressures were referenced to the sphere pressure. 

The rocket-motor combustion- chamber pressure was measured by a transducer with 
a range of 1000 pounds per square Inch tapped into the end of the combustion 
chamber. 

Thermocouples and calorimeters .- All thermocouples used in this investiga- 
tion were chromel-alumel. Two wire sizes, 0.005 inch and 0.010 inch in diam- 
eter, were chosen as a compromise between anticipated response and durability 
requirements. The larger wire thermocouples were generally used in locations 
close to the nozzle exit. 

Drawings showing the manner in which the thermocouple probes and calorim- 
eters were constructed are presented as figure 4. The chromel and alumel wires 
for the thermocouple probes were joined by first flattening the end of the two 
wires slightly and then spot-welding together. The junction thus formed was 
then dressed with a file to roughly the same diameter as the wire. The ther- 
mocouple wires were then passed through a 2-inch length of two-hole ceramic 


5 


tubing and spliced to larger gage chromel and alumel wires. This assembly was 
then slipped into a 3-inch length of 3/l6-inch outside diameter stainless steel 
tubing and sealed at both ends with ceramic cement. Construction of the calo- 
rimeters was similar with the exception that the formation of the thermocouple 
junction and its contact with the bottom surface of the copper disk was accom- 
plished with one spot-welding operation. The weight of each copper disk was 
determined by means of a chemical balance before attachment of the thermocouple 

wires. Tables listing the property parameter — 22- and 22 f or the individ- 

12 4 A 

ual thermocouples and calorimeters and their locations on the reflection plate 
for the model configurations investigated are presented in figure 5* Also 
shown in this figure are enlarged silhouettes of two typical thermocouple 
junctions. 


Data Reduction and Accuracy 

All data obtained during a test run were recorded with recording oscillo- 
graphs. These records were then transcribed manually into punch-card form and 
machine processed. Data readout was in the form of absolute and differential 
pressures and temperatures. 

Experimental heat flux to the disk calorimeters and thermocouples was cal- 
culated from the temperature-time histories by using the following equation: 

g _ 

A ~ ^At 

d wc Wc t 

where k = "jj - for the thermocouples and k = ~ for the calorimeters. (For 

the thermocouples, A was taken as the entire surface area of the exposed ther- 
mocouple wire. ) The rate of change of temperature with time was obtained manu- 
ally from the temperature-time plots; this procedure introduces a random error 
in the computed heating rates. The specific heat of the thermocouples was 
taken as the arithmetical average between the specific heats for the chromel 
and alumel wires. 

The use of equation (l) involves the following assumptions: (l) heat loss 

due to radiation is negligible; (2) heat loss due to conduction through the 
thermocouple lead wires is negligible; ( 3 ) the thermal conductivity of the disk 
calorimeters and themocouples is such that the temperature difference through 
the calorimeter disk and through the exposed thermocouple wire is negligible; 
and (4) the heat-transfer coefficient is invariant with temperature. For the 
present investigation the error in the measured heat flux introduced by the 
first two assumptions was calculated to total about 6 percent (less than 2 per- 
cent by radiation and about ^ percent by conduction) for the time period of 
1.2 seconds after motor ignition for which heat-transfer coefficients are pre- 
sented. These losses were calculated for the thermocouple having the highest 
temperature rise in this time period. This error was considered to be small 
for this investigation; therefore, the heat-flux calculations were not cor- 
rected for these losses. 
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Another possible source of error in the temperature data is the heat flux 
due to radiation from the exhaust plume to the thermocouples. This error was 
found to be negligible from a preliminary test run in which the exhaust plume 
barely impinged upon the reflection plate (essentially no backflow). The tem- 
perature rise experienced by the thermocouples for this run was essentially the 
same as the temperature rise of the sphere thermocouple which could not see the 
exhaust plume. 

Heat-transfer coefficients were obtained from a modified definition of the 
convective heat-transfer coefficient 


h = 


ci/A 

T t - T 


( 2 ) 


where the usual recovery or adiabatic wall temperature has been replaced by a 
total temperature. The adiabatic wall temperature was impractical to evaluate 
in the present investigation and, since it is not readily available in most 
cases, the heat-transf er coefficient was based on the total temperature of the 
backflow. This total temperature was, in turn, assumed to be equal to the 
rocket-motor combustion- chamber temperature (4800° F) . 

The accuracy of the temperature and pressure measurements is estimated to 
be as follows: 


Temperature, °F ±3 

Combustion- chamber pressure, Ib/sq in ±10 

Static pressure, lb/sq in ±0.001 

Sphere pressure, lb/sq in ±0.0004 


For the reflection-plate static pressures and the sphere pressure, the error 
listed is the maximum estimated to occur when only the accuracy of the pressure 
transducer and the accuracy in reading the oscillograph traces are considered. 
The total-pressure probe readings were subject to an unknown error because of 
possible misalinement with the flow. The probe was alined parallel with the 
reflection-plate center line for all configurations. 


RESUITS AND DISCUSSION 


Pressure Data 

Sphere pressure and rocket-motor performance .- The sphere pressure, rocket- 
motor combustion- chamber pressure, and rocket-motor thrust measured during the 
several rocket-motor firings of this investigation are presented in figures 6 
and 7- The geometric altitude corresponding to the sphere pressure scale in 
figure 6 is the right-hand ordinate in this figure. The useful burning time of 
the rocket motor was about k seconds except in the case of the motor of configu- 
ration II (fig. 7)* For this configuration the useful burning time was about 
3-5 seconds as indicated by the rapid decrease in combustion- chamber pressure 
at this time value. The curves in figure 7 "between t = 0 and t - 0.2 have 
been arbitrarily faired through the data scatter obtained in this time period. 
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It will be noted in figures 6 and 7 that there is a lack of consistency 
between the sets of curves. For example, the rocket motor with the highest 
combustion- chamber pressure did not have the highest value of thrust and the 
motor with the highest thrust did not yield the greater sphere pressure. It 
is felt that the main reasons for this inconsistency are the probable differ- 
ences in throat areas and combustion temperatures of the propellants between 
the several nozzles used in this investigation. A graphite insert was used 
for the throat of the nozzle and there was no strict dimensional control over 
these inserts. 


Nozzle pressure ratio and exhaust boundary.- The variation of nozzle total- 
pressure ratio with time presented in figure 8 was computed by using values of 
p^ c and p^ midway between the maximum and minimum values at a given value of 

time t in figures 6 and 7* The approximate nozzle static-pressure ratio 
(ratio of nozzle exit pressure to the sphere pressure) for the same conditions 


may be obtained from figure 8 by using the nominal value 


plotted values of 


^t,c 

Poo ' 


Pe 

*t,c 


0 . 005 and the 


The angle (3 at which the computed exhaust boundary strikes the reflection 
plate for the two spacing ratios is given in figure 8 for two nozzle pressure 
ratios. The exhaust boundaries were computed by the method of reference 7 by 
using the median nozzle pressure ratio curve of figure 8, a ratio of specific 
heats of 1.22, and the nozzle dimensions. 


Reflection-plate static pr essures.- Static pressures measured on the reflec- 
tion plate are presented in the form of a ratio to the sphere pressure for sev- 
eral periods of time after motor ignition in figures 9 a nd 10. Figure 9 pre- 
sents this ratio for a line perpendicular to the vertical center line (bend) of 
the reflection plate and located 8 inches ahead or upstream of the nozzle exit. 
Similarly, figure 10 shows the static pressure distributions obtained along a 
line parallel to the bend in the reflection plate. This line was offset about 
1 inch from the bend. The curves of figures 9 and 10 were obtained by first 
plotting the absolute pressure time history for each orifice and fairing a 
smooth curve through these data points. The ratio p/p^ was then computed 

from these faired curves and the corresponding sphere pressure data. The pres- 
sure distributions were then obtained from these plots of p/p against time 

' 00 

for the indicated times in figures 9 and 10 by assuming flow symmetry about the 
plate center line; that is, pressure orifice Pp was transposed to a location 

between p^ and p^. The pressure data obtained for t < 1.0 second suffers 
from lack of accuracy because of the low pressures generally obtained in this 
time period. In some cases negative absolute pressures were indicated in this 
time period. For these cases the curve for pressure plotted against time was 
faired to zero pressure at t = 0. 

If, in figure 9> only that portion of the pres sure -distribution curves 
between the plate center line and 8 inches from the center line is considered 
for t > 1.0, a region of high pressure is noted between 5 and 7 inches from 
the plate center line for both spacing ratios. The pressures measured at the 
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outer edge of the plate for the larger spacing ratio (fig. 9(*0) appear to he 
abnormally high when compared with the other measured plate pressures. The 
pressure distributions obtained parallel to the plate center line for 

— = 8.02 (fig. 10(a)) show a rapid rise in pressure from the exit plane to 
d e 

the line of orifices of figure 9 (8 inches upstream from the nozzle exit plane) 
for t > 0.4. However , this rise in pressure occurs only at t = 0.4 for the 
larger spacing configuration of figure 10(b) . For other values of t, the pres- 
sure increases gradually as the distance increases. Also to be noted in fig- 
ures 9 and 10 is the increase in V B /v m for a given orifice with time. Fur- 
ther discussion of these observations will be made in the section entitled 
"Flow Model. " 


Backflow total pressures .- The total pressure of the exhaust backflow was 
measured at three locations on the reflection plate for the smaller nozzle 
spacing configuration and at one location for the larger spacing configuration. 
(See fig. 3(b).) Total -pres sure time histories for these locations are pre- 
sented in figure 11 as a ratio of the sphere pressure. Also presented in this 
figure are ref lection -plate static-pressure-ratio time histories corresponding 
to the total-pressure time histories. Static pressures shown for configura- 
tions I and II were obtained for the probe location from reflection-plate 
static -pres sure contours. For configurations III and IV, however , it was nec- 
essary to use the static pressures from an orifice located some distance from 
the probe (see figs. ll(c) and 11(d)); thus, the pressure at this orifice may 
or may not be Indicative of the static pressure at the probe location. It will 
be noted that the probe total -pres sure ratio eventually becomes less than the 
plate static-pressure ratio for all configurations except configuration III. 
This difference would Indicate that the backflow has ceased and the probe Is 
subject to flow from the opposite direction. 


Backf low Mach numbers and unit Reynolds number .- The interpretation of 
measured total pressures in a rarefied gas is complicated by viscous forces 
which become appreciable compared with the inertia or pressure forces in the 
flow; that is, the total pressure is related to a Reynolds number in addition 
to the Mach number and static pressure. A discussion of the analytical work 
of reference 8 in reference 9 shows that, for the total-pressure probe consid- 
ered in reference 9 j the difference between a measured total pressure and the 
total pressure determined from the Rayleigh pitot formula for known values of 
Mach number and static pressure is large at low supersonic Mach numbers. In 
reference 9, the measured value exceeded the computed value by 60 percent at 
Mach 1.5* At hypersonic Mach numbers the difference approaches a near constant 
value of 11 percent. 


The analysis of reference 8 was used in the calculation of the backflow 
Mach number from the pressure data of figure 11. The resulting Mach number 
variation with time for these data is presented in figure 12. It will be noted 
here that there Is a considerable variation in Mach number with time as well as 

laterally across the plate for a given time, (configurations I, II, and III 
are for ^ = 8.02.j 
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In order to compute the unit Reynolds number of the backflow, it was nec- 
essary to evaluate the viscosity of the exhaust gas. The composition of the 
exhaust gas mixture and its gas constant R were known from data supplied by 
the propellant manufacturer. The viscosities of the individual gases were 
obtained from reference 10 and the viscosity of the mixture was then approxi- 
mated by combining the percent by weight of the viscosity of the individual 
gases. The use of a more accurate method of determining the viscosity of the 
mixture, such as the procedure of reference 11, involves further assumptions 
that did not appear to be justified from the data available. 

The results of the Reynolds number calculations are presented in the lower 
portion of figure 12. These Reynolds number time histories correspond to the 
Mach number time histories shown in the upper portion of figure 12. The unit 
Reynolds number is seen to be of the order of magnitude of 1 so that a Reynolds 
number based on the diameter of the thermocouple wire is extremely small. Cal- 
culation of the Knudsen number K from these values of Mach number and Reynolds 
number by the relationship 


K = 1.26/7 -M- 

u Re,T 

yielded a minimum value of K of about 38. Thus, the data obtained for the 
thermocouples were well into the free-molecule flow region. (See ref. 13.) 


Temperature Data 

Sphere gas temperature . - Temperature time histories obtained for the sphere 
thermocouple for the five configurations of this investigation are presented in 
figure 15* This thermocouple was located behind the reflection plate and could 
not "see" the rocket motor or the exhaust plume. The temperature rise indi- 
cated by the sphere thermocouple in figure 13 for a given configuration results 
from the hot exhaust gases and the compression of these gases within the sphere 
during the burning time of the rocket motor. Any of the plate thermocouples 
not subject to the exhaust backflow should also register this temperature rise. 
It will be noted that the sphere gas temperature rise for configurations II 
and V was somewhat lower than the temperature rise for the other three 
configurations . 

Position T1 thermocouple .- The temperature time histories for the thermo- 
couple In position T1 (see fig. 3(b)) for the five configurations are pre- 
sented in figure l4. These time histories are characterized by an initial, 
rapid rise in temperature followed by a leveling off or slight decrease in tem- 
perature. The smaller nozzle spacing configurations (configurations I, II, 
and III) exhibit the greatest overall temperature rise, this rise being two to 

three times as great as that for the — = 17-05 configurations (configura- 
te 

tions IV and V) . Since this thermocouple was located In the plane of the noz- 
zle exit, it was subject to direct exhaust-gas impingement (especially so for 
the smaller nozzle spacing) and greater radiation effects from the exhaust 
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gases. Direct exhaust-gas impingement during the early moments after motor 
ignition would account for the rapid temperature rise noted in this time 
period. The decrease in the rate of temperature rise or the leveling off of 
the temperature between 1.6 and 2.4 seconds, depending upon the configuration, 
is thought to be due to changing flow conditions as the sphere pressure 
increases. That is, the thermocouple is subject initially to direct impinge- 
ment but as the sphere pressure increases, the plume boundary strikes the 
reflection plate further downstream and the thermocouple comes under the influ- 
ence of the backflow. The differences in thermocouple response between the 
configurations for a given nozzle spacing can be attributed to the differences 
in height of the thermocouple above the reflection plate. The zero thermo- 
couple height of configurations II and III may place the thermocouple in closer 
proximity to any high-pressure region on the plate because of plume impingement 
than the 1-irich height of configuration I. For the larger nozzle spacing, how- 
ever, by assuming little or no direct impingement, the zero-height thermocouple 
of configuration V would be subject to backflow that has been cooled by passage 
over the reflection plate. 

In view of the uncertainty as to the flow conditions existing for thermo- 
couple Tl, these data will not be discussed further insofar as backflow heat- 
transfer coefficients are concerned. However, as a matter of interest, the 
order to magnitude of h for configurations II and III is about 1 X 10-3 
Btu/ft2-sec-°F. 

Plate thermocouples and disk calorimeters .- Temperature time histories of 
the thermocouples and disk calorimeters located in positions T2 to Tl6 are pre- 
sented in figures 15 to 19 for the five configurations. The temperature time 
history of the thermocouple used to measure the sphere gas temperature is also 
presented in each figure for comparison purposes. 


Of general interest in figures 15 to 19 is the response of the plate 
thermocouples relative to that for the sphere thermocouple. It will be noted 

4- j-*. 4 - ti- T ( ^t __ Q pk)i ,, J _ -L-U _ 4- -I ^,4 _ 


that only for configuration I 


= 1 inch | is the temperature rise 


of the plate thermocouples generally greater than that for the sphere thermo- 
couple over the burning time of the rocket motor. For configurations II and III 
(y T = 0) at this nozzle spacing, the temperature rise of the plate thermocouples 

is generally less than that for the sphere thermocouple except near the exit 
plane of the nozzle. This reduction in temperature rise for these configura- 
tions is attributed to cooling effects of the reflection plate on the backflow 
next to the plate. The distance perpendicular to the plate that this cooling 
effect (or thermal gradient) would extend would increase as the distance 
upstream from the nozzle exit increased. Therefore, it is reasonable to expect 
that even those thermocouples that were extended 1 inch above the plate surface 
(configurations I and IV) and located some distance upstream of the nozzle may 
have been adversely influenced by the presence of the plate. At the larger 
nozzle spacing the plate thermocouples exhibit an initial temperature rise that 
is generally greater than that for the sphere thermocouple. However, a rather 
abrupt change in the temperature rise occurs between about 1.5 to 2.0 seconds 
after motor ignition for most of the thermocouples and the temperature of these 
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thermocouples thereafter increases at a much slower rate or shows a decrease in 
temperature. This temperature response of the plate thermocouples suggests a 
significant change in the "backflow pattern or complete cessation of the backflow 
rather than attainment of equilibrium temperature by the thermocouples. 

The temperature rise of the disk calorimeters noted in figures 15 to 19 is, 
of course, less than that for the thermocouples because of the greater heat 
capacity of the calorimeters. In figures 17(a) and 18(a) it will be noted that 
a calorimeter was mounted perpendicular to the backflow at position T10 for con- 
figuration III and at position T5 for configuration IV. 

Validity of nozzl e cluster simulation method.- The significance of the 
pressure and temperature data discussed in the foregoing paragraphs is dependent 
upon whether the reflection plate was a valid representation of the planes of 
interference that would normally occur for a cluster of three nozzles. Data and 
schlieren photographs from reference 12 and an unpublished source yielded strong 
evidence that the reflection plate would be a valid representation of these 
interference planes. Part of this evidence is shown in the schlieren photo- 
graphs of figure 20. The top photograph in this figure is from reference 12 and 
depicts a sonic nozzle operating at a static-pressure ratio of about 5 X 10^. 

The bottom photograph is a schlieren view of a cluster of four sonic nozzles 
operating at a static pressure ratio of about 20. It will be noted in this fig- 
ure that the shock produced at the impingement point of the exhaust flows from 
two adjacent nozzles (bottom photograph) is very similar to the one produced by 
the flow from a single nozzle impinging upon a flat plate (top photograph) . 
Base-pressure measurements on the four-nozzle cluster in figure 20 show that 
the presence of a reflection plate to simulate the flow from two of the nozzles 
had no significant effect on base pressures for nozzle spaeings greater than 
C + 

— = 0-75- At lesser spaeings the boundary layer on the plate caused an effec- 
de 

tive nozzle spacing less than the geometric spacing. 

It would appear from the foregoing that the use of a reflection plate to 
simulate a nozzle cluster was a valid experimental procedure for the purpose 
of this investigation. However, it should be pointed out that it is possible 
to conceive of circumstances where erroneous data may result from use of this 
method. For example, it can be postulated that in a space environment the mean 
free path of the molecules in the outermost portions of the exhaust plumes may 
be so large that the plumes would initially intermingle without mutual interfer- 
ence. In this case, the backflow would not originate along the lines of initial 
plume contact but from a higher pressure region within the merged plumes where 
formation of interference shocks would be possible. For this situation the 
approaching streamlines from the individual plumes would meet at a shallower 
angle, and thus reduce the quantity of backflow or, at the extreme, eliminate it 
entirely. Therefore, in this circumstance the reflection plate would produce 
a greater amount of backflow than would occur from the actual cluster of nozzles 
being simulated or, if viewed from the standpoint of nozzle spacing, the reflec- 
tion plate would create an effective spacing ratio smaller than the geometric 
one. Since it is not known whether this or any other similar circumstance 
existed for the present investigation, the data are presented on the basis that 
the reflection plate was a valid experimental procedure for the purpose intended. 


12 


Flow Model 


The supersonic velocities measured in the backflow of this investigation 
support the results of reference 6 wherein data are presented that indicate that 
the backflow is supersonic a short distance from the common impingement point of 
the nozzle plumes. Therefore, for the balance of this discussion it is assumed 
that the backflow becomes supersonic at or very near the line of impingement of 
the exhaust plume on the reflection plate. With this in mind and consideration 
of the data presented and the post-fire photograph of the reflection plate in 
figure 2l(a) , a possible explanation of the behavior of the backflow in the 
present investigation can be set forth. A factor which enters into the explana- 
tion is the circulation of the gases within the sphere induced by the rocket 
exhaust flow. This circulating flow and the other factors under discussion are 
shown in the drawing of figure 2 l(b) . In this figure the supporting framework 
for the reflection plate is not shown and only the top portion of the rocket 
motor case is depicted. 

In figure 2l(a) a sharply defined semicircular area on the reflection plate 
will be noted wherein the dust from the rocket grain deposited on the plate is 
less than that outside this area. The region in which the dust "ring" crosses 
the line of static-pressure orifices coincides with the high-pressure area noted 
in figure 9 ( a ). This condition suggests that the dust ring (also shown in 
fig. 21(b)) represents a region of high pressure relative to the sphere pressure 
and that the increase in pressure is due to a sudden decrease in velocity of the 
backflow. The latter deduction follows from the lack of dust deposit in a 
region of known supersonic flow and the heavy deposit in the area surrounding 
this region. The question arises here as to whether the dust ring represents 
the interface between the induced flow and the backflow or whether it represents 
a deceleration of the backflow to subsonic speeds by a shock such as would occur 
for overexpanded flow in a nozzle and the interface occurs beyond the dust ring 
as shown in figure 2 l(b) . The latter supposition is felt to be the case here 
as two of the three calorimeters (t8 and T 9 ) located just beyond the dust ring 
show heating rates comparable to those applied to the thermocouples within the 
area enclosed by the dust ring. (See fig. 22 .) Thermocouples located further 
beyond the dust ring on the lower portion of the plate, positions Til to Tl6, 
show an overall temperature increase that is not much greater or somewhat less 
than that for the sphere thermocouple (fig. 16(b)). The heating rates for these 
thermocouples shown in figure 22 are larger than the rate for the sphere thermo- 
couple initially but show a sudden decrease to values about equal to or less 
than the sphere thermocouple rate soon after motor ignition. This decrease sug- 
gests that the induced flow becomes predominant in this region of the plate as 
shown In figure 21 (b) soon after motor ignition and the less than sphere thermo- 
couple heating rates result from cooling of the induced flow as it passes over 
the sphere wall and floor before reaching the model. The rather abrupt decrease 
in the temperature rise of the thermocouples of the larger nozzle spacing con- 
figurations (figs. 18 and 19) suggests that the induced flow becomes a dominant 
fact over much of the reflection plate soon after motor ignition because of the 
reduced amount of backflow for this spacing ratio. The concept of an interface 
between the induced flow and backflow leads to the deduction that the interface 
acts as an effective base which moves closer to the exit as the sphere pressure 
increases. If this were the case, it could account for the increase in 
reflection-plate pressures with time noted previously for figures 9 and 10. 
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Heat-Transfer Coefficients 


Heat-transfer- coefficie nt profiles . - Average heat-transfer-coefficient dis- 
tributions for small-diameter cylinders (thermocouples) placed in the backflow 
for the five configurations are shown pictorially in figure 23- These coeffi- 
cients are an arithmetical average of the values obtained at five time intervals 
from t = 0.^ to t = 1.2 seconds. Only the data to 1.2 seconds were used as 
these data were believed to be free of any induced-flow effects. A time period 
of 0.4 second after motor ignition was allowed for establishment of the flow. 
Average heat-transfer coefficients obtained for the disk calorimeters are listed 
in the upper left corner of figure 23. It will be noted in this figure that 
only half of the reflection plate is depicted. The line of thermocouples on the 
other half (line emanating at 22^° with respect to the plate center line from 

position T1 in fig. 3(a)) lias been transposed to the half shown. Although the 
temperature data for the thermocouple at position T1 were not considered further 
because of possible direct exhaust-plume impingement, the curves for the distri- 
butions have been faired with the knowledge that, in any case, the value of h 
would be greater at position T1 than at any other location. 

The effect of the thermocouple height above the reflection plate on the 
heat-transfer coefficient is evident in figure 23 "by comparing configuration. I 
with configurations II and III and configuration IV with V. For configuration V 
the heat meters were at yq = 1.0 inch, whereas the thermocouples were at 

yrp = 0. Comparison of configuration II with configuration III gives an indica- 
tion of the repeatability of the data since these configurations differed only 
in the locations of the calorimeters and the total -pres sure probe. The heat 
meters mounted perpendicular to the backflow (figs. 23(c) and 23(d)) show 
expected higher values of average heat-transfer coefficients. For all config- 
urations the coefficients decrease with increasing distance from the exit, the 
coefficients closest to the exit having a maximum value of about 3 X 10”^ 

Btu/sq. ft-sec-°R. 

Heat-transfer-coefficient correlation . - The average heat-transfer coeffi- 
cients obtained for the thermocouples at positions T4, T5, and T 6 are compared 
with empirical and theoretical values for cylinders from references 13 and l4 
in figure 24. In this figure the heat-transfer coefficient is presented as a 
Nusselt number variation with Reynolds number with the viscosity and thermal 
conductivity in these parameters evaluated at stagnation conditions and the den- 
sity evaluated at free-stream conditions. The values of h used in the calcu- 
lation of % u ^ 0 were taken from the profiles of figure 23* It was assumed for 

the model = 8.02 that the velocities obtained for configurations II and III 
d e 

were applicable to configuration I in order to compute NRe,o for the thermo- 
couple positions listed in this figure. Since the values of h used were based 
on the backflow total temperature, the values of Nifu^o obtained were somewhat 
higher than would be the case if h were based on the recovery temperature; 
that is, in free-molecule flow the recovery or adiabatic wall temperature is 



greater than the total temperature. (For the molecular speed ratios M of 
the present Investigation, the recovery factor is about 1.3* ) 

It Is seen in figure 24 that there is disagreement between the experimental 
values and the work of reference 13 . (The curves of this reference have been 
extended by about an order of magnitude smaller in fig. 2k.) For agreement, the 
experimental values need to be about two orders of magnitude larger in Np e ,o 
or about one order of magnitude smaller in N^ u 0 . Extrapolation of the data of 
reference l4 to the Np e ,o values of the present Investigation shows good 
agreement. However, in view of the amount and the range of the data of refer- 
ence 13 , this agreement appears to be coincidental. It is felt that the primary 
cause for the disagreement in the data of figure 2k is the value of NR e ,o- 
Uncertainty exists as to the correctness of the calculated values of this param- 
eter because of the number of assumptions and corrections Involved. As sug- 
gested in reference 13> the static pressures as well as the total pressures may 
be subject to errors due to viscous effects. In any case, however, large 
changes in the static pressure or Mach number would be required to increase 
%te,o "k° magnitude required for agreement with the data of reference 13- 


CONCLUDING REMARKS 


Data obtained from an exploratory investigation of exhaust backflow 
resulting from the mutual interference of exhaust gases from a simulated cluster 
of three nozzles operating in a near-space environment indicated that backflow 
was present with spacings of 8.02 and 17*05 nozzle exit diameters between noz- 
zles. Pressure data obtained in the backflow region near the nozzle exit showed 
that the backflow was supersonic and that the flow over small-diameter thermo- 
couples placed in the backflow region was well into the free-molecule flow 
region. Average heat-transfer coefficients obtained for small-diameter cylin- 
ders (thermocouples) in the backflow region were greater than those predicted by 
empirical and theoretical methods for cylinders In high-speed rarefied gas flow. 
This result may have been due, in part, to the inability to calculate an accu- 
rate Reynolds number from the pressure data available for the backflow region. 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Administration, 
Va., May 27, 1965. 
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Reflection plate 



Figure L- Sketch showing manner in which a reflection plate was used to simulate the interference from exhaust boundaries of adjacent rocket 

motors. Figure not drawn to scale. 
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(a) Sketches showing location and typical installation of instrumentation. 
Figure 3.- Reflection plate. 
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Figure 4.- Sketch showing thermocouple and calorimeter construction. Figure not drawn to scale. 




Figure 5.- Enlarged silhouettes of typical 
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Figure 10.- Concluded. 
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(b) Configuration II. 
Figure 23 .- Continued. 




(c) Configuration III. 
Figure 23.- Continued. 
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(d) Configuration IV. 
Figure 23.- Continued. 





(e) Configuration V. 
Figure 23.- Concluded. 
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